Doped AB 2 Hubbard Chain: Spiral, Nagaoka and RVB States, Phase Separation and 

Luttinger Liquid Behavior 
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We present an extensive numerical study of the Hubbard model on the doped AB2 chain, both in 
the weak coupling and the infinite-U limit. Due to the special unit cell topology, this system displays 
a rich variety of phases as function of hole doping (8) away from half-filling. Near half-filling, spiral 
states develop in the weak coupling regime, while Nagaoka itinerant ferromagnetism is observed in 
the infinite-U limit. For higher doping the system phase-separates before reaching a Mott insulating 
phase of short-range RVB states at S = 1/3. Moreover, for S > 1/3 we observe a crossover, which 
anticipates the Luttinger liquid behavior for 8 > 2/3. 

PACS numbers: 71.10.Fd, 74.20.Mn, 75.30.Kz 
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Low-dimensional strongly correlated electron systems 
have attracted great attention in the last two decades. 
The reason dates back to Anderson's proposal that the 
t — J version of the Hubbard model might carry the basic 
mechanisms underlying the high-Tc superconductivity 
observed in CuC>2 compounds. Despite that this remains 
an open issue, the above suggestion fertilized intensive in- 
vestigations on many related fundamental topics, such us 
itinerant electron magnetism, Mott metal-insulator tran- 
sitions and quantum critical phenomena. Amongst sev- 
eral features of interest, we mention the possibility of re- 
alization of spiral*, Nagaoka-*** and resonating-valcncc- 
bond (RVB) states^, spatially separated phases 7 ** and 
Luttinger liquid behavior—, which may present strong de- 
viations from the Landau Fermi liquid theory. 

In this work, we report numerical results of the Hub- 
bard model on the doped AB2 chain away from half 
filling, which show that its special unit cell topology 
greatly enriches the phase diagram found in the doped 
standard linear chain. In fact, all features mentioned 
above are shown to be associated with well defined 
ground state (GS) phases of this doped chain. Doped 
AB2-Hubbard chains were previously studied through 
Hartree-Fock, quantum monte carlo and exact diagonal- 
ization (ED) techniques both in the weak and strong 
coupling Unfits**, including also the t — J model** us- 
ing density matrix renormalization group (DMRG) and 
recurrent variational Ansatzes, and the infinite-U limit** 
using ED. In particular, these chains represent an al- 
ternative route to reaching two-dimensional quantum 
physics from one-dimensional system o 11 ' 13 . At half fill- 
ing the AB2-Hubbard chain exhibits a quantum ferri- 
magnetic GS 10 ' 14 i 15 ' 17 , whose magnetic excitations have 
been studied in detail both in the weak and strong cou- 
pling limits* 7 -, and in the light of the quantum Heisen- 
berg mode l 17 ' 18 . Further studies have considered the 
anisotropic** and isotropic*** critical behavior of the 
AB2-quantum-Heisenberg model, including its spheri- 
cal version**, and the statistical mechanics of the AB2- 



On the experimental side, the AB2 chain topology is 
of relevance to the understanding of the physics of some 
low-dimensional strongly correlated electronic systems. 
One class is the line of trimer clusters present in fos- 
fates with formula A 3 Cu3(P04) 4! where A=Ca***&**i**, 
gr 24 1 25 i 26 i 27 and Pb***&**. The trimers have three Cu 2+ 
(S = 1/2) paramagnetic ions antiferromagnetically cou- 
pled. Although the superexchange intcrtrimer interac- 
tion is much weaker than the intratrimcr coupling, it 
proves sufficient to turn them bulk fcrrimagncts. Another 
quasi-one-dimcnsional inorganic material closely associ- 
ated with the ferrimagnetic phase of the AB2 chain is the 
NiCu bimetallic chain**. These compounds display alter- 
nating Ni 2+ (S = 1) and Cu 2+ (S = 1/2) ions connected 
through suitable ligands in a line; and are modeled by 
the alternating spin- ^ /spin- 1 antifcrromagnetic Heisen- 
berg chain**. We also would like to mention a more re- 
cently synthesized organic ferrimagnetic compound con- 
sisting of three S — 1/2 paramagnetic radicals*- 1 in its 
magnetic unit cell, as well as possible connections with 
the physics of the oxocuprates**. 



This paper is organized as follows: in Sec. II we in- 
troduce the model system and the numerical techniques 
used to calculate several quantities suitable to charac- 
terize the occurrence of distinct phases as function of 
doping and Coulomb coupling. In Sec. Ill, we discuss 
spiral and Nagaoka states at low hole doping, whose mag- 
netic properties are shown to exhibit very interesting fea- 
tures in the weak and infinite-U limit, respectively. In 
Sec. IV we show that for higher hole doping the system 
phase separates, before reaching a Mott insulating phase 
of short-range RVB states at 6 = 1/3. In Sec. V we dis- 
cuss several features of the crossover region, which takes 
place before the Luttinger liquid behavior observed for 
5 > 2/3. Finally, in Sec. VI we present a summary and 
some conclusions concerning the reported results. 



II. MODEL DESCRIPTION AND METHODS 

The AB2 chain is a bipartite lattice with three sites 
(named A, Hi and B 2 ) per unit cell, as illustrated in Fig. 
QJa). The Hubbard Hamiltonian for a lattice with N c 
unit cells and N sites reads: 

N a N 

H = -tV2 [&L(^ + ^+i, CT ) + H.c.] + [/^n lT n 4 , 

l=l,o i=l 

(1) 

where A\ a and b\ a = -^(B\ la + B\ M ) are the creation 
operators of an electron with spin a at site A and in 
a bonding state between sites B\ and B2 of the cell I, 
respectively, t(= 1) is the hopping amplitude and U is 
the Coulomb coupling. For U = 00, double occupancy is 
completely excluded and the Hamiltonian takes the form: 

H = -tV2 PoHriAia + A l+1 .„) + H.C.]P G , (2) 

i=l,<7 

where Pg = ~~ n i] n i[) IS the Gutzwiller projector 

operator. The model is invariant under the interchange 
of the B sites of the same cell, a symmetry that implies 
in a well defined local parity (pj = ±1) for the GS wave 
function. As a result, in computing some quantities we 
find it convenient to use the effective linear chain (ELC) 
generated by the map illustrated in Figs. Ufa) andQTb), 
i. e., any quantity Xb,i associated with a B site at cell I 
of the ELC is given by Xb x ,i + Xb 2 ,1- This mapping does 
not change the physical content of the GS and excited 
states, being used only to expose in a more clear fashion 
some properties of these states. 

In the tight-binding description (U = 0) this model 
presents three bands^: one flat with vY c odd parity states 
[antibonding orbitals, a\ a = -j={b\ 1(J — B\ la )] and en- 
ergy e = 0; and two dispersive branches, 

e±(&) = ±2\/2cos(fc/2), (3) 

with k = 2nl/N Cl I = 0, 1, 2, N c - 1, built from A 
sites and bonding (even parity) orbitals, as shown in Fig. 
[Tic). At half filling (N e = N, where N e is the number 
of electrons) the GS total spin S g is degenerate, with 
S g ranging from the minimum value (0 or 1/2) to S g — 
\Nb — Na\/2, where Na(Nb) is the number of sites in 
the A (B) sublattice. As proved by Lieb^ the Coulomb 
repulsion lifts this huge degeneracy and selects the 

S g = \N B -N A \/2 = S Lleb (4) 

ground state for any finite U, giving rise to a ferrimag- 
netic GS^^l 

On the other hand, for U = 00, one hole (N e = N — 1) 
and periodic boundary conditions (BC's), the system sat- 
isfies the requirements of Nagaoka's theorem for satu- 
rated fcrromagnctis m 10 ! 12 . For Nagaoka ferromagnetism 
and Lieb ferrimagnetism the GS is homogeneous in par- 
ity with pi = — 1 for any cell I. Due to this symmetry, 
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FIG. 1: (Color online), (a) Illustration of the AB2 chain 
showing A, B\ and B2 sites; (b) Illustration of the effective 
linear chain (ELC). (c) Electronic bands of the tight-binding 
model: two dispersive (continuous line) and one flat (dashed 
line) . 



the spectrum of the AB2 chain in the Heisenberg limit 
(U >> t, N e = N) at the sector p = —1 is identical to 
that of the alternating Heisenberg spin-i/spin-1 chain^S. 

Here we focus on the effect of hole doping, 5 = 
1 — (N e /N), both in the weak coupling and the infinite- 
U limit, using exact diagonalization (ED) through the 
Lanczos algorithm for closed BC's and DMRG for open 
BC's^. In the ED procedure, the BC's are such to mini- 
mize the energy, except for U = 2 and S < 1/3 [Fig. [UJc)] 
in which the BC's (periodic or antiperiodic) are such that 
the Fermi wave vector fcp in the thermodynamic limit is 
included in the set of wave vectors for the finite system^. 
We used finite size DMRG for open chains with A sites 
in its extrcma, keeping 364 to 546 states per block in 
the last sweep. The maximum discarded weight in the 
last sweep was typically ~ 10 -7 , except for odd phases 
and U = 2, where the discarded weight was ~ 10~ 5 . In 
the DMRG calculations we treated B\ and B2 as a com- 
posite site with 9 states for U = 00 and 16 states for 
U = 2. However, by considering the parity symmetry, 
we can decompose this supersite into the two possible 
symmetry sectors +1 and —1 . Within this scheme, we 
have considered all parity symmetry sectors of the form 
( — ) X {+) N "~ X , with x contiguous cells of odd parity in 
one side of the open chain and N c — x contiguous cells 
of even parity in the other. In addition, we have verified 
the stability of this phase separation against the forma- 
tion of a mixed phase composed of smaller domains. The 
energy is studied as function of x for increasing number 
of states kept per block in order to localize the value of x 
for which the energy is minimum, as shown in Figs. [2Ja) 
and[Uh). The phase-separated boundaries are thus de- 
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termined by the limiting dopings for which an inhomoge- 
neous phase (non-uniform parities) is observed. We have 
also developed a simple variational approach for U = oo 
and 5 < 1/3, which is explained in detail in Appendix A. 
The results calculated using this approach are shown in 
Figs. Hid) and G^c) . 

In Fig. Ogc) (U = 2) and Fig. C^d) (U = oo) we present 
the average parity, 



(5) 



as function of doping, computed using the above- 
mentioned methods. In both regimes, we observe the 
occurrence of an homogeneous phase near half filling with 
p = —1. For higher doping, i. e., <5pg(f7) < S < (1/3) 
[<5pg(2) fts 0.07 and <5pg(oo) w 0.22] the system phase 
separates in one region with odd parity cells and the other 
with even ones. For S > 1/3 the GS is homogeneous with 
p = l. 

In order to present an overview of the conducting prop- 
erties of the AB2 chain phases in the infinite- U limit, we 
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FIG. 3: (a) and (b): Magnetic structure factor for U — 2 
and N = 100 in the underdoped region, (c) Total spin per 
cell S/N c as function of 5 for U = 00. 



display in Fig. fJIe) the quantity^ 

, 2irqi 



(6) 



2N r 



calculated in the ELC using ED, where L 
jrij, rij is the electron density at site j and q is such that 
^ = 2, with p and q co-primes. The phase of z^ cor- 
responds to the GS expectation value of the position op- 
erator, while its modulus defines the localization length; 
in an insulator, \z^\ — » 1, as L — ► 00, while in a conduc- 
tor, \z^\ — > 0, for closed boundary conditions^. The 
increase of \z^\ with system size for 5 = 2/3 and 1/3, 
as well as in the phase separated region, are evidences 
of insulating phases at these dopings. These conclusions 
will be better fundamented by studying the Drude weight 
using ED and the charge gap for larger systems with 
DMRG. 



III. SPIRAL STATES AND SATURATED 
FERROMAGNETISM 

In Figs. [3ja) and[3]Jb) we display the magnetic struc- 
ture factor 



FIG. 2: DMRG results for the energy difference between the 
lowest energy at the symmetry sector ( — ) x (+) Na ~ x and the 
GS energy for N = 100 (N c = 33) for (a) U = 2 at <5 = 0.32 
(triangle up), S = 0.26 (diamond) and 5 — 0.18 (triangle 
down); (b) U = 00 at <5 = 0.32 (triangle up), S = 0.28 (trian- 
gle up) and S = 0.24 (triangle down), taking 108 states per 
block (open symbols) and 216 states per block (filled sym- 
bols). Average local parity p as function of S for (c) U = 2 
and (d) U = oo. (e) ED results for \z^\. Dashed lines are 
guides to the eye. 



2Nc+l 



S(q) 



?Lieb{Sh 



-TT E ^- m) <S*-S m ), (7) 

Lieb + 1) , 



calculated at S z = and U = 2 using DMRG for the 
ELC. First, notice the presence of peaks at q = and 
q = 7r revealing the ferrimagnetic order at half filling. 
These peaks sustain up to two holes (S = 0.02); however, 
it is not clear whether the ferrimagnetic phase is robust 
against doping in the thermodynamic limit. Indeed, by 
increasing the hole doping, spiral peaks at ^-dependent 
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FIG. 4: (Color online). ED results of (a) N c times the energy 
gap Ao between the saturated ferromagnetism energy ($ = 
$f) and the lowest energy state for an Aharanov Bohn flux $ 
as function of <£> = |$ — $f| at 8 = 1/6 and N c = 4 (dashed- 
dot line), 6 (dashed line) and 8 (solid line). ED calculation for 
the <3> dependent behavior of the (b) Spin correlation function 
(S ce ]i (Zo) ■ S ce ji (io + 1)) between the cell spins as function of I 
and (c) magnetic structure factor as function of lattice wave 
vector q at 8 — 1/6. (d) Charge structure factor calculated at 
the lowest energy state for any $ at 8 = 1/6 for N c = 4 (•), 
6 (A) and 8 (V). 



positions appear near q = and q = tt. The analysis of 
the charge gap, 



E{N e + 1) + E(N e - 1) - 2E(N e ), 



(8) 



suggests that these states are metallic, in opposition to 
the Mott insulating ferrimagnctic state at 8 = 0. It is 
worth mentioning that the occurrence of spiral phases in 
oxocuprates has been a challenging and topical subject 3 ^. 

In Fig. [3^c) we present the GS total spin as function 
of doping for U = oo. For S < <5pg(oo) itinerant satu- 
rated ferromagnetism due to hole kinematics (Nagaoka 
mechanism) is observed. It is interesting to notice that 
our estimate for the upper hole density (> 0.2) beyond 
which Nagaoka ferromagnetism is unstable is in very 
good agreement with similar predictions for ladder a 37 ! 38 
and the square lattice 3 ^. 

We have also considered the presence of an Aharonov- 
Bohn flux $ for a closed chain through the gauge trans- 
formation: 



ho 



Aire*** 1 '*: 



(9) 



with $o = hc/e = 1. The flux variation is equivalent 
to a change in the boundary condition: $ = represents 
periodic and $ = 1/2 antipcriodic boundary conditions. 
In Fig. H^a) we present the dependence of the energy 
gap Ao between the lowest energy state for a flux $ and 
that for saturated ferromagnetism ($ = as function 
of $ = |$ - $f| at 5 = 1/6. We have identified many 



level crossings in this curve. In fact, as the flux increases 
from <f>p, the total spin decreases from the maximum 
value, S = N e /2, to the minimum value S — (S = 1/2) 
for N e even (odd) , a behavior also observed in the square 
lattice^. Notice that N c Aq tends to saturation with sys- 
tem size, indicating that the level spacings decrease with 
1/N C . These results suggest that the thermodynamic GS 
displays spontaneously SU(2) symmetry breaking as a re- 
sult of an ergodic combination of infinitely many states 
(N c — > oo), including the singlet spiral stated. In Figs. 
[4fb) and Htc) we present the spin correlation function 
between cell spins S ce u(l) = Sa(1) +Sb 1 (Z) + Sb 2 {1) and 
the magnetic structure factor 



S(q) 



J_ J2 e lq{l " m) (Sod, (I) ■ Scdi (m)) (10) 
C (l,m) 



as function of distance I and wave vector q = 2nl/N Cl 
I = 0, ...,N C , respectively. As we can observe, the satu- 
rated ferromagnetic and the spiral singlet states are adi- 
abatically connected, such that all states contributing to 
the thermodynamic GS exhibit long-range ordering. In 
particular, as the flux increases from <&f the peak of S(q) 
at q = (saturated ferromagnetism) steadily decreases, 
while the spiral state peak at q = 2tt/N c increases. We 
noted also that the charge structure factor 



(11) 



<Z,m) 



where An; = ni — (ni) and ni is the electron occupation 
number at cell I, is not affected by the flux variation 
and displays a peak at 2kp = tt [Fig. HJd)], where fcp 
is the tight-binding spinless Fermi wave vector—, with 
k F = 3tt(5, 5 < 1/3. 



IV. PHASE SEPARATION AND RVB STATES 

In the phase-separated regime the charge compressibil- 
ity diverges following the linear dependence of the energy 
with doping. In Figs. [5]and[n]we present some properties 
of the GS in this regime calculated through DMRG for 
the ELC. First we notice that all these properties clearly 
exhibit some modulation on the same sublattice in the 
metallic odd parity region due to charge itinerancy. In 
particular, this modulation is stronger in the U = 2 spiral 
phase as evidenced by the correlation function (Si • S») 
shown in Fig. [5^a), but also noticed in the itinerant Na- 
gaoka phase (U = oo) as manifested by the site magne- 
tization (Sf) shown in Fig. [5^c). On the other hand, 
in the insulating even parity phase a flat behavior is ob- 
served, except for boundary and interface effects. These 
paramagnetic phases [see Figs. E^b) and [5fcl)] are char- 
acterized by strong singlet correlations between spins at 
sites Bi and B 2 at the same cell, i. e., (Sb 1 Sb 2 ) ~ —0.20 
(pa —0.41) for U = 2 (= oo), as shown in Figs. [5jb) and 
(5jd). In contrast, in the metallic phase this correlation 



■5 



varies very little with U and indicates robust triplet cor- 
relations, i. e., (S Bl • S B . 2 ) w 0.13 (w 0.16) for U = 2 
(= oo). Notice that in the absence of hole hopping, even 
when restricted to a cell as in the insulating phase, the 
value of (Sbj • Sb 2 ) in a singlet (triplet) state should be 
—0.75 (0.25). The hole density (nh,i) is shown in Figs, 
[(^a) and^b). In the odd parity metallic phase, holes do 
not occupy antibonding orbitals, whereas in the even par- 
ity insulating phase these orbitals are accessible for them. 
Therefore, in the first case the hole densities at sites A 
and Bi + B 2 are very similar. This may also occur in the 
second case if double occupancy is excluded (U = oo). 
An illustration of the phase-separated regime for U = oo 
is shown in Fig. [6f c) . In this coupling limit, unsaturated 
ferromagnetism was suggested to occur in ladders 3 ^ and 
the square lattice 3 ^ as an intermediate phase between sat- 
urated ferromagnetism and paramagnetism as function of 
doping. However, in the context of the t — J model the 
situation is more complex and predictions of phase sep- 
aration, both for ladders 3 ^ and the square lattice ^ 42 ' 43 , 
and stripe formation for the square lattice^ 2 - have been 
reported. 

At 6 = 1/3, i. e., one hole per A site for open BC's 
using DMRGii, the GS has even parity and is fully 
dominated by the Mott insulating phase (even parity) 
illustrated in Fig. [5{c) for U = oo. The charge gap 
A c = wh ere (J.+ = [E(N e +AN e )-E(N e )]/AN e , 

AN e > (AN e /N -> 0), and ^_ = E(N e ) - E(N e - 1), 
must be calculated with care. First, notice that adding 
electrons to S = 1/3 places the system in the phase- 
separated (inhomogeneous) region where the chemical 
potential /i is flat. Indeed, by comparing results using 
DMRG and ED calculations for U = oo, for which A c 
presents little finite size corrections [Fig. HJa)], we con- 
cluded that boundary effects are minimized by taking 
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FIG. 5: GS properties at 5 = 0.18 (U = 2) and S = 0.28 
(U = oo) for N = 100 using DMRG. (a) Spin correlation 
function (Si ■ Si) for U = 2. (b) Expectation value of Sf 
for U = oo in the sector S z — S g . Spin correlation function 
(Sbj ■ Ss 2 )i for (c) U = 2 and (d) U = oo. — (+) indicates 
odd (even) local parity. Dashed lines are guides to the eye. 
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FIG. 6: (Color online). GS properties at 5 = 0.18 (U = 2) and 
8 = 0.28 {U = oo) for N = 100 using DMRG. Expectation 
value of rih,i for (a) U = 2 and (b) U — oo. Effective linear 
chain notation:(») identifies A sites and (o) B\ + Bi at the 
same cell, (c) Illustration of the GS for U — oo in the phase- 
separated regime: singlet bonds are represented by ellipses 
and holes by circles. — (+) indicates odd (even) local parity. 
Dashed lines are guides to the eye. 



AN e = 2 and placing the symmetry inverted cells at the 
chain center. We thus find [Fig. 0[a)] A c w 0.21 (w 0.96) 
for U = 2 (= oo). This problem is absent in the case of 
hole doping since the phase is homogeneous. The extrap- 
olated spin gap, 



A s = E(S = 1) - E(S = 0), 



(12) 



characterized by symmetry inversion of a cell at the chain 
center, is also shown in Fig. [T^a) for U = 2 (As ~ 0.18) 
and U = oo (As ~ 0.16), with the spin gap at U = oo 
presenting little finite size dependence. It is a quite mas- 
sive excitation with the magnon localized at the odd sym- 
metry cell, mostly at the B sites, as shown in Fig. [TJb). 
In this context, Sierra et aL— found As « 0.27 using the 
t - J model (J = U 2 /U) for J = 0.35t, i. e., U w 11.43. 
We have confirmed this result by studying the U depen- 
dence of As using ED. In Fig. [5] we show that the 
spin correlation functions at 5 = 1/3, calculated using 
DMRG, present a fast decay and can be fitted with the 
exponential form exp [—{I — l c )/£)], where £ is the corre- 
lation length, I is the cell index in the ELC and l c denotes 
the central cell of the system. This behavior is expected 
from the presence of a finite spin gap. The values of 
i for the correlations (S A (l c ) ■ S A (l)}, (S A (l c ) • S B (Z)) 
and (S B (lc) ■ S B (0) are w 0.4 (2.2), 0.25 (0.45) and 0.39 
(0.75), respectively, for U = oo (U — 2), with l c de- 
noting the central cell. Thus, except for the correla- 
tion (S A (l c ) ■ S A (l)) at U = 2, the correlation length 
is extremely short with spins correlated only within a 
cell. Further, the calculated bulk values of (Sb 1 ■ Sb 2 ) 
at S = 1/3 are in very good agreement with those in the 
even phase of the separated region shown in Figs. [SJb) 
and [5]^d) . The above results support a short-range- RVB 
(SR-RVB)ii state for the GS at 5 = 1/3, as illustrated in 
Fig. [5Jd). In this context, Sierra et al*^ reached similar 
conclusions using the t— J model on the AB2 chain, while 
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FIG. 7: DMRG results for the (a) size dependence of the 
charge (A c ) and spin (As) gaps as function of 1/N C at 
S — 1/3: solid lines are polynomial fittings. ED results for 
the (b) expectation values of S\(l), S%(1) and the correlation 
function (Ss 1 (Z) ■ Ss 2 (Z)) at spin sector S z = 1 as function of 
cell number I. The ± signs below the horizontal axis in (b) 
indicate the cell parity. 



Giesekus has proved^ that a SR-RVB state is the GS of 
a non-bipartite lattice with the same local symmetry but 
a different hopping pattern. 
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FIG. 8: Spin correlation functions (a) (Sa(Z c ) ■ Sa(/)), ( b ) 
(S A (/ C ) • Ss(0) and (c) (S B (Z C ) • Su(Z)) as function of Z - Z c 
in the ELC at 5 = 1/3 for JV C = 33 using DMRG; in the 
above expressions l c denotes the central cell; (0) refers to 
(7 = 2 and (♦) to U = do. Dashed lines are guides to the 
eye. (d) Illustration of the GS at S = 1/3, singlet bonds are 
represented by ellipses and holes by circles. 



V. LUTTINGER LIQUID BEHAVIOR 

We now focus on the behavior of the system for 1/3 < 
5 < 1 by considering a chain with closed boundary con- 
ditions and N c = 8 for U = oo using ED. The first no- 
ticeable feature is the behavior of the spin correlation 
functions after doping the (5=1/3 GS with two holes. 
The value of (Sb 1 (Iq) ■ Sb 2 (Iq)) (where Iq denotes an 
arbitrary cell) changes from -0.41 to -0.28. This varia- 
tion can be understood by considering that the two holes 
added to the system break two singlet bonds and reside 
predominately at B sites. In this picture the correla- 



tion function would amounts to 



N 



N, 



K-0.41) 



-0.31, 



which is close to -0.28. Furthermore, the spin correlation 
functions shown in Fig. 02a) evidence the formation of 
long ranged bonds between electrons on B sites, while 
the other correlations remain short ranged, as in the 
5=1/3 ground state. This fact indicates that the elec- 
trons picked from the SR-RVB by hole doping arc antifer- 
romagnetically coupled and delocalized through the sys- 
tem, as illustrated in Fig. EJb). In order to describe the 
system behavior for finite dopings, we display in Fig. [He) 
the correlation function (Sb ± (Jo) • Sb 2 (Iq)) and electronic 
densities as function of 5. Notice that for 1/3 < 5 < 2/3 
the electronic density at A sites is almost fixed, while 
that at B sites are monotonically depopulated. As a 
consequence, (Sb^o) ■ Sb 2 Go)) continuously vanishes as 
the doping increases. Moreover, in Fig. [H^d) we show 
the relevant nearest-neighbor spin correlation functions. 
These correlations display quite different magnitudes at 



S = 1/3, but their values approach each other for S > 2/3. 
We thus consider the doping interval l/3<<5<2/3as 
a crossover region, where doping starts to build the Lut- 
tinger liquid which is fully established for 5 > 2/3. 

We have also calculated the charge compressibility k 
through 

X = 4- = T^ N " + 2 ) + E ( N < - 2 ) - 2£ W]> ( 13 ) 
UqK 4 

where V is the volume and no = ^f- is the electronic 
density; the charge excitation velocity 



EjAk, S = 0) - E GS 
Ak 



(14) 



with Afc = 2ir/L and L the system length; and the Drude 
weight 



D 



L 

47T 



d 2 E($) 



(15) 



where <& m in is the flux value that minimizes the energy^. 
In an insulating phase these quantities satisfy the limits 
below 



lim 



N c 



X 
u p 
D 



oo; 
oo; 
0. 



(16) 
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FIG. 9: (a) ED results for the indicated spin correlation func- 
tions as function of cell index I, with S_b = Ssj + S_b 2 . (b) 
Illustration of the GS at S — 1/3 doped with two holes: 
singlet bonds are represented by ellipses and holes by cir- 
cles, (c) ED results for the spin correlation functions be- 
tween B sites at the same cell (Ss 1 (Zo) ■ Ss 2 (Zo)), electron 
densities at A sites (tm.(Zo)) and at B = B\ + B2 sites 
(rifl(Zo)) = (ws 1 (Zo) + tib 2 (Io))- (d) ED results for the indi- 
cated nearest-neighbor spin correlation functions as function 
of S, with Sb = + S_b 2 - In ( a ), ( c ) an d (d) Zo denotes an 
arbitrary cell. 



while for a metal x, u p and D are finite. As shown in 
Fig. [lOl at S = 1/3, x an d u p increases, while D decreases 
with system size for both U = 2 and U — 00, although 
the insulating character is better evidenced for U = 00 
due to its sizable charge gap, as shown in Fig. [7]Ja). At 
the other commensurate density, 5 = 2/3, we can see the 
signals of an insulating phase for U = 00, while for U = 2 
we does not observe any especial behavior. In order to 
clarify this point, we have used DMRG to study the size 
dependence of the charge gap for larger systems at this 
doping. For a finite open chain, the occupation of two 
holes per cell tends to 5 = 2/3 in the thermodynamic 
limit. In Fig. [TlTa). we can clearly observe that for 
U = 00 the system is in a Mott insulating phase with 
A c w 0.15; however, the gap for U = 2 is extremely 
small. In order to better understand the U-dependcncc 
of this gap, we have also calculated A c for intermediate 
values of U, as also shown in Fig. [Tlja). In the inset of 
Fig. [TTT a) we have fit A C (U) using an expression similar 
to the limiting behavior of the charge gap as U — > of 
the Lieb-Wu solution for a linear chain at half filling 47 : 
U a exp (—b/x), in which a 0.61 and b ss 7.95 are fitting 




FIG. 10: ED results for the charge susceptibility x, the charge 
excitation velocity u p and the Drude weight D, for [7 = 2 
[(a),(b) and (c)] and U = 00 [(d), (e) and (f)], and iV c = 4 
(•), 6 (A) and 8 (T). 



parameters. Notice, however, that contrary to the Lieb- 
Wu solution^ 7 -, A c saturates to a finite value (w 0.15) for 
U = 00. On the other hand, similarly to the linear chain 
at half filling^, the data shown in Fig. [TlTb) indicates 
the absence of spin gap at S = 2/3 in the thermodynamic 
limit for both [7 = 2 and U = 00. 

In the Luttinger model, it is well known^ that \, u p 



and D arc related through 



with 



D 



Kn 



2u p K p , 



TtU c 



2x 



(17) 



(18) 



where K p is the exponent governing the decay of the 
correlation functions. In order to probe the doped region 
for which the lower energy spectrum of the AB2 chain 
can be mapped onto the Luttinger model, we consider 
the ratio 



Ratio 



v^xA' 



(19) 



which must be equal to one if the system is in the LL 
universality class«. 

Since the AB2 chain is not strictly one-dimensional, 
care must be taken with the length scales (V and L) 
in Eqs. $13j) , dHJ) and (HJ). For U = 0, the orbitals at 
sites A and bonding orbitals at sites B are translationally 
equivalent and both build the dispersive branches shown 
in Fig. [Tic). In this case, the system can be mapped 
onto a tight-binding linear chain with 2Af c sites and a 
rescaled hopping parameter, t — ► t\f 7 2 1 with K p = 1. In 
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FIG. 11: DMRG results for the (a) Charge gap A c as function 
of 1/N C at 5 = 2/3 using DMRG; the inset presents extrapo- 
lated values of the charge gap as function of U . DMRG cal- 
culation of the (b) Spin gap Ag as function of 1/N C for U — 2 
(•) and U = oo (♦). Solid lines are polynomial fittings, except 
in the inset of (a), where we have used an essential singularity 
form as explained in the text. 



order that Eq. (fT5| matches this result for ep < 0, we 
must choose V = L = 2N C with e(fc) = — 2\/2cos(fc); 
or, likewise, V = L = N c and the dispersions as written 
in Eq. ([3]). In both cases kp = ^n^, with n$ = 
Consider, for example, the former option. For U = 0, the 
charge excitation velocity is equal to the Fermi velocity 
up, which can be easily calculated as 



Uf 



de(k) 



dk 



= 2\/2sm(fc F ). 



k — A; j 



On the other hand, substituting the GS energy, 

-8V2- 



Ecs{n ) 



-7V c sin(-n ), 



(20) 



(21) 



into the continuous version of Eq. (|13[) . we obtain, 
1 d 2 E GS 



X 



V 



dnl 



7r\/2 sin(— no) 



(22) 
(23) 

TBI) we find, as 



Using now Eqs. ([20]) and (J23J) in Eq, 
expected, K p = 1. 

We now turn to the interacting case using ED. As 
shown in Fig. [TWal the LL character is quite clear 
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FIG. 12: (a) ED results for the Ra,tio=u p /s/Dx/n. (b) ED 



results for K p as function of 8. 




FIG. 13: Spin correlation functions C(l) for (a) U = 2 and 
(b) U = 00 at <5 = 88/106 for N = 106 using DMRG: solid 
lines are fittings using Eq. (f25|) . 



for 5 > 2/3, while for 1/3 < 5 < 2/3 we identify the 
crossover region. The ED results for K p are presented in 
Fig- ITST b) . Notice that K p is close to 1 (non- interacting 
fermions) for U = 2; while, K p is close to 1/2 (non- 
interacting spinless fermions) for U = oo^d^. In order to 
check these results, we used DMRG to calculate the ELC 
spin correlation function 



C(l) 



s \i-j\,i 



(24) 



whose asymptotic behavior should match that for the 
Luttinger model^: 



C LL (l) 



cos(2M)[ln(0] 1/2 



(25) 



In Eq. ([24]) we have considered an average over all pos- 
sible pairs of sites separated by the same distance I, a 
procedure that reduces open boundary effects. In Figs. 
□Ha) andCCJb) we show C(l) calculated at 5 = 88/106 
for U = 2 and U = 00, respectively. Also shown are 
the fittings to C(l) using Cll{1) with kp = ^n and K p 
taken from the results shown in Fig. [T27 b) after linear in- 



terpolation: K p 



0. 



(U = 2) and K p = 0.57 (U = 00). 



Motivated by a compromise between large values of I and 
minimum boundary effects, we have considered interme- 
diate values of I in the fitting, which is quite good for 
both values of U. We thus conclude that the Luttinger 
model correctly describes the low energy physics of the 
AB 2 chain for 8 > 2/3. 
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VI. SUMMARY AND CONCLUSIONS 

In summary, the numerical results presented here have 
clearly evidenced the rich phase diagram exhibited by 
the Hubbard model on the doped AB2 chain both for 
U = 2 and in the infinite- U limit. We have shown that 
at the commensurate dopings (5=1/3 and 2/3 the sys- 
tem display insulating phases, although for U = 2 the 
charge gap A c is very small at 5 = 2/3, with indica- 
tions that A c present an essential singularity as U —* 0. 
For U = 2 and 5 < 0.02 the GS exhibit a ferrimag- 
netic phase reminiscent of the undoped regime, while for 
0.02 < 5 < 0.07 incommensurate magnetic correlations 
are observed. For U = 00 and 5 = the GS total spin 
is degenerate, whereas for < 5 < 0.225 hole itinerancy 
(Nagaoka mechanism) sets a fully polarized GS. In this 
case, we have also observed the presence of an extensive 
number of low-lying levels with total spin ranging from 
the minimum value to S rnax — 1 and level spacing decaying 
with system size as 1/N C . For higher doping, the system 
phase separates into coexisting metallic and insulating 
phases for Sps{U) < 6 < 1/3 (with Sps{oo) w 0.225 and 
<5ps(2) w 0.07). The insulating state presents a finite 
spin gap and fully fills the system at 5 = 1/3, which is 
well described by a short-ranged-RVB state. Finally, a 
crossover region is observed for 1/3 < 5 < 2/3, while a 
Luttinger liquid behavior is explicitly characterized for 
5 > 2/3. 

In closing, we would like to stress that the above- 
reported results might also stimulate further exper- 
imental and theoretical investigations on quasi-one- 
dimensional compounds displaying complex unit cell 
structures^. 

We acknowledge useful discussions with A. L. Malvezzi 
and M. H. Oliveira. This work was supported by CNPq, 
Finep, FACEPE and CAPES (Brazilian agencies). 



APPENDIX A: VARIATIONAL APPROACH FOR 

U = 00 AND S < 1/3 

In the metallic saturated ferromagnetic region (parity 
symmetry -1) the energy as function of doping is known 




0,2 0,4 0,6 0,: 
X 



FIG. 14: Energy per unit cell e as function of the fraction 
x = L-/N c for S = 0.28. The physical region is also shown. 



to have a non-interacting spinless fermion behavior 
4V2 



-L- sin(fc F _/2), 



(Al) 



where kp_ = nVh—, fh— = Nh-/L—, and L_ is the lin- 
ear size of the system. On the other hand, in the insulat- 
ing paramagnetic phase (SR-RVB states with even parity 
symmetry) at 5 = 1/3 (one hole per cell) 



Nh 



L 



+ • 



(A2) 



and the energy per cell e+ is almost independent of the 
system linear size and can be estimated either by using 
ED or DMRG: 



-2.021. 



(A3) 



Let us now consider a phase-separated regime in which 
a paramagnetic phase with size L + coexists with a ferro- 
magnetic one with size £_ , so the energy per cell reads 



Nr. 



4V2L_ . ,tt 



It is convenient to write i>h- as 

N h -{N C -L-) 35 



Vh- = 



(1-x) 



L- 



(A4) 



(A5) 



where N h = N h+ + N h -, N c = L + + L-, x = L-/N c 
and N = 3N C . Using the above notation we rewrite Eq. 
(I All) in the form below 



„ x 4^2 . 

e(x) = (1 — x)e + ism 



3(5- 1 



(A6) 



Here we should notice the presence of a singularity at 
x = for any finite value of 5 ^ 1/3 (sec Fig. [T4| . 
However, the region of physical values of x is defined by 



< N h+ < N h , i. e., 



1 - 35 < x < 1. 



(A7) 



(A8) 



In Fig. [14] we present e(x) for 8 = 0.28, in which the 
physical region is 0.16 < x < 1 and can be found by Eq. 
(IA8|) . with a minimum in e(x) for x 0.49. 

The value of x which minimizes the energy for a given 



5, x = x(5), satisfies the equation 
can be written as 



dc(x) 



dx 



4V2 



cos(y) + ysin(y), 



where 



it 35- 1 



0, which 



(A9) 



(A10) 
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The roots of Eq. 
conduct to 



9j are numerically calculated and 



x = 1, for 
x« 3.071 -9.213<5, for 



5 < 0.225; 
0.225 < 5< |. 



(All) 



We thus conclude that 5ps(oo) « 0.225, which is in very 
good agreement with ED and DMRG calculations. 

The magnetization is null at the even phase and max- 
imum at the odd one. We can thus derive the following 
expression for the GS total spin per unit cell: 



Nr. 



1 

2N t 
1 
2 



{N e - 2L+) (A12) 
[3(1 — 8) — 2(1 - x)] (A13) 



The dependence of the average parity p on 6 can also be 
easily written as 



p = 1 - 2x 



(A14) 



Finally, using Eq. (jAlip for x, the above results for p 
and S g arc plotted in Figs. [2](d) and [He), respectively, 
and shown to be in excellent agreement with the ED and 
DMRG calculations. 
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